Phase diagram of the metal-insulator transition in 2D electronic systems 
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We investigated the interdependence of the effects of disorder and carrier correlations on the 
metal-insulator transition in two-dimensional electronic systems. We present a quantitative metal- 
insulator phase diagram. Depending on the carrier density we find two different types of metal- 
insulator transition - a continuous localization for r 3 < 8 and a discontinuous transition at higher 
r s . The critical level of disorder at the transition decreases with decreasing carrier density. At very 
low carrier densities we find that the system is always insulating. The value of the conductivity 
at the transition is consistent with recent experimental measurements. The self-consistent method 
which we have developed includes the effects of both disorder and correlations on the transition, 
using a density relaxation theory with the Coulomb correlations determined from numerical simu- 
lation data. 

PACS numbers: 73.20.Dx,7f .30.+h,73.40.-c 
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A metal-insulator transition (MIT) has now been ob- 
served in a number of different 2D electronic systems 
over a wide range of carrier densities and levels of dis- 
order. At the transition the carrier density parameter 
r s , which measures the strength of the carrier correla- 
tions, covers values 7 < r s < 35. The strength of 
disorder at the transition also covers a wide range, and 
depends on the value of the critical r s . In Si/SiGe || the 
transition occurs at low carrier densities and high mobil- 
ities. For the transition in Si MOSFETs which occurs at 
higher densities, the mobilities are much smaller. A re- 
cent experimental phase diagram Q provides a relation 
between the strength of the correlations and the disor- 
der. It shows that the critical level of the disorder at the 
transition diminishes as the correlations grow stronger. 

In spite of a great deal of experimental and theoreti- 
cal work the nature of the conducting state remains un- 
clear and controversial. Unlike conventional metals, its 
Coulomb interaction energy is typically an order of mag- 
nitude larger than its Fermi energy. Numerous proposals 
have been made about the causes of the stabilization of 
this conducting phase. These include strong Coulomb 
repulsions between the carriers |||q], or an anomalous 
enhancement of the spin-orbit interaction due to the bro- 
ken inversion symmetry of the confining potential well in 
Si MOSFETs @. References g] proposed that the con- 
ducting phase is superconducting, with the pairing of the 
carriers being mediated by the dynamic correlation hole 
surrounding each carrier. Thakur and Neilson |l demon- 
strated the existence of superconducting pairing due to 
strong Coulomb correlations in the presence of disorder. 
They found up to levels of disorder typical of high quality 
Si MOSFETs that the superconductivity persists. 

The conducting phase is destroyed not only at low 
but also at high carrier densities. Kravchenko et al |l]] 
originally observed the insulator-to-metal transition by 
increasing the carrier density. Recently it has been re- 
ported that if the carrier density is increased further there 
is an upper critical value where the conducting phase 



again becomes unstable and the system enters another 
insulating phase [Q. This insulating phase has the char- 
acteristics of a single-particle localized state, while the 
insulating phase at low densities is a coherent insulator 
with properties similar to a Wigner crystal or glass. The 
critical density for the second transition, typically around 
r s < 7, is still not low enough for electron interactions to 
be neglected. In this paper we treat Coulomb interactions 
and disorder on an equal footing in order to investigate 
the transitions at both small and large r s . 

In the strong correlation limit and in the absence of 
disorder, electrons localize to form a Wigner solid with 
long range crystalline order. Here we consider interact- 
ing charge carriers in the presence of weak disorder which 
would destroy any long range order. We have previously 
proposed that strong Coulomb correlations in the pres- 
ence of disorder can localize electrons into a coherent 
glassy insulator with liquid-like short-range order || . We 
studied the metal-insulator transition in terms of the er- 
godic to nonergodic transition of the charge density fluc- 
tuations. The possibility of localization into such a glassy 
state has also been discussed recently by Chakravarty et 
al @. 

We define the order parameter for the glass state as 
f(q) =linii_ >00 *(g,t), where $(q,t) = (N(q,t)\N(q,0)) 
is the Kubo-relaxation function for the dynamical den- 
sity variable N(q,t). The normalized variable N(q,t) = 
p(q,t) I \/x(q)i where p(q,t) is the usual density fluctua- 
tion operator and xil) is the static susceptibility. When 
the order parameter is non-zero, spontaneous density 
fluctuations do not decay at infinite time and the sys- 
tem will be an insulator. Conversely if f(q) is zero then 
our system is in a conducting phase. Since our order pa- 
rameter would be zero for any conducting phase it gives 
no indication about the precise nature of the conducting 
phase. 

Within the Mori-Zwanzig formalism jll| $(q,t) is cal- 
culated in terms of the memory function M(q, t), which 
we evaluate using mode-coupling theory []12[ . In the limit 
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t — > oo the relaxation function reduces to 
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m = T 



n(q)/M(q) 



(1) 



where fl(q) = q 2 /(m*x(q)) and M(q) = limt^,*, M(q, t). 
m* is the carrier effective mass. We express M (q) = 
M cc (q) + M ic (q), where 

M cc (q) = 2^-2 E ' + V(\q - q'\) 

<;' 
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M ic (q) = K<|^ inip (g)| 2 > + <|W surf (g)| 2 >] 
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The M cc (q) part of the memory function originates from 
the interactions between the carriers, and the Mi C (q) 
from the carrier-disorder interactions. V(q) — 2ire 2 /eq is 
the Coulomb potential with dielectric constant e. Ui mp (q) 
is the carrier-impurity potential for randomly distributed 
monovalent Coulombic impurities of density n.j which are 
embedded in the carrier layer. W SUI t(q) is the surface 
roughness scattering term. Details of the disorder poten- 
tials used are given in Ref. M . 

To evaluate the static susceptibility x(l) with correla- 
tions, we use the generalized Random Phase Approxima- 
tion expression, 



X(q) = *>(«)/{! + V(q)(l - G(q)) X o(q)h 



(3) 



where Xo{q) is the Lindhard function for non- interacting 
electrons. The local field factor G(q) accounts for the 
correlations between the carriers. We evaluate G(q) from 
ground state properties of the electron liquid jl3[ using 
the fluctuation-dissipation theorem |Q . 

The memory functions M cc (q) and Mj c (g) mutually 
influence each other through f(q). The nature of the lo- 
calized state at the transition is largely determined by 
which of these memory functions is dominant. At low 
densities and small levels of disorder the interactions be- 
tween the carriers dominate and M cc (q) is much larger 
than Mi c (q). In this case the localization is primarily 
caused by many-body effects and the localized state is a 
coherent frozen insulator. At the transition the order pa- 
rameter f(q) jumps from zero to non-zero and the system 
undergoes a discontinuous transition. In contrast, at high 
densities and high levels of disorder, where the carrier- 
disorder scattering dominates, Mi C {q) is much larger than 
M cc (q). The localization transition in this case is to a 
non-coherent state where the carriers localize indepen- 
dently. We find that this transition is continuous, with 
the order parameter f(q) continuously increasing with 
disorder. Our model thus predicts the two distinct types 
of metal-insulator transition. 






FIG. 1. Order parameter f(q). Curve labels are impurity 
densities rii. Surface roughness is included in all cases. 

a. r a = 8. rii is in units of 10 10 cm -2 . 

b. r s = 12. rii is in units of 10 9 cm -2 . /(g) is zero for 
in < 1.92 x 10 9 cm~ 2 . 

c. r s — 18. rii is in units of 10 8 cm -2 . Dashed line is rii = 
with no surface roughness. 
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We solve Eqs. [l] and || self-consistently for the order 
parameter f(q). We vary the strength of the disorder 
potential by varying the impurity density rii while keep- 
ing the surface roughness at each r s fixed. 

In Fig. |l| we show the order parameter f(q) for three 
values of r s = 8, 12 and 18. The multiple curves for 
each fixed r s correspond to increasing levels of disorder. 
At the highest carrier density, Fig. |j]a, the localization is 
in general dominated by the Mi C (q) part of the memory 
function. Here f{q) increases continuously with impurity 
density, indicating continuous localization. These char- 
acteristic features that f(q) evolves continuously with the 
disorder level and is non-zero for any disorder indicates a 
localized state where the particles localize independently 
This is similar to Anderson-type localization. Since f(q) 
is always non-zero there is no conducting phase when 
r s = 8. The peak in f(q) at q = 2kp is a result of the 
well-known cusp in the two-dimensional xo(q) and is of 
no importance here. 

In contrast, in Fig. for r s — 12 the order parameter 
is zero for small non-zero levels of disorder. This indi- 
cates the existence of a conducting phase. If we increase 
the level of disorder, a critical value is passed at which 
f(q) jumps discontinuously to non-zero values and at that 
point there is a metal-insulator transition || . The short- 
range coherent order of the insulator state is reflected 
by a peak in f(q) at the reciprocal nearest neighbor dis- 
tance, q k, 2Akp. As we increase the disorder further, 
the f{q) increases continuously. For very high levels of 
disorder, rii > 6 x 10 9 cm~ 2 , the overall shape of the 
f(q) evolves towards a Gaussian indicating the develop- 
ment of a non-coherent insulator. We find in the range 
8 < r s < 18 that the critical level of disorder required 
to localize the carriers decreases with increasing r s and 
that the non-coherent state also occurs at decreasing rii. 

By r s — 18 the system localizes without disorder and 
there is no conducting phase. This large r s localization 
is drivenpurely by the correlations between the carriers. 
In Fig. |ljc f(q) is non-zero for rii — and no surface 
roughness scattering (dashed line) . Without disorder the 
f(q) goes to zero as q goes to zero. The solid lines show 
f(q) with surface roughness scattering included. 

The property that the critical level of disorder for the 
metal-insulator transition is dependent on the carrier 
density is associated with the changing strength of the 
correlations. The zero temperature phase diagram in Fig. 
H shows the relationship between the carrier density and 
the critical level of disorder at the transition. Disorder 
includes both the impurity density m and surface rough- 
ness. For r s < 8 the system is insulating in the presence 
of any disorder, and the order parameter f(q) is always 
non-zero. For 8 < r s < 18 the order parameter is zero for 
rii less than a critical value and we have the conducting 
phase. At the right phase boundary the transition to the 
insulating phase is discontinuous and f(q) jumps discon- 
tinuously to non-zero values as we cross it. The critical 
level of disorder for the transition decreases with decreas- 
ing density. Near the phase boundary where the disorder 



level is slightly larger than the critical value the system 
is in the coherent insulating state. For very large disor- 
der the f(q) has evolved into a Gaussian-like shape (see 
Fig. 0b) and the system becomes a non-coherent insula- 
tor. For r s > 18 the order parameter is non-zero without 
disorder and the conducting phase has disappeared. Our 
phase diagram has common features with the conceptual 
phase diagram discussed by Chakravarty et al [ fL0| . 
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FIG. 2. Zero temperature phase diagram. Axes are im- 
purity density m and carrier density parameter r s . In the 
insulating phases the order parameter f(q) > 0. In the con- 
ducting phase f(q) = 0. The transition at higher densities 
is continuous. When r s > 9 the f(q) discontinuously jumps 
from zero to non-zero values at the transition. Inset shows 
the corresponding scattering rate ~//(2Ef) at the transition 
as a function of ru. 

The inset in Fig. ^| shows the relation between rii and 
the scattering rate 7 off the disorder. Within the memory 
function formalism the non-linear equation for 7 is given 

by© 

i7=-^^E« 2 M^imp(?)| 2 } + <|WW<7)| 2 }] 

x ( x(g) V Mi,h) r41 
W°Hq)J i + i7<M<M7)/x (0) (9) ' 1 1 

where ^0(5,17) = (V^) [x (0) (<7>i7) -X (0) (<7)] is the re- 
laxation spectrum for non-interacting carriers that scat- 
ter off the disorder. Thakur and Neilson || showed for 
the density range 5 < r s < 10 that the superconduct- 
ing phase persists at least up to values of r y/(2Ep) = 1. 
They also showed that the effective potential becomes 
more strongly attractive as r s increases. This suggests 
that our conducting phase in Fig. |^, which is only found 
for ^/{2Ef) % 1, is a superconductor. 

Using the Drude model, 7 can be related to the zero 
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temperature conductivity, a = (ne 2 / 'm*)7~ 1 . In Fig. || 
we show a at the transition as a function of r s . In the 
range 9 < r s < 14, a lies between 1 < a/(e 2 /h) < 3. 
This is consistent with experimental values. At lower 
densities the a calculated within the Drude model in- 
creases rapidly. This is associated with the decrease in 
the critical level of disorder at the transition. Since the 
Drude approximation is inapplicable for large r s this in- 
crease in a at low density should not necessarily be seen 
exp er iment ally. 
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FIG. 3. Conductivity a at the critical disorder for the 
transition as a function of r s . 
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